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NORTH ATLANTIC HURRICANE CHANGES 



Abstract. Baycsian statistical models were developed for the 
number of tropical cyclones and the rate at which these cyclones 
became hurricanes in the North Atlantic. We find that, controlling 
for the cold tongue index and the North Atlantic oscillation index, 
there is high probability that the number of cyclones has increased 
in the past thirty years; but the rate at which these storms become 
hurricanes appears to be constant. We also investigate storm inten- 
sity by measuring the distribution of individual storm lifetime in 
days, storm track length, and Emanuel's power dissiptation index. 
We find little evidence that the distribution of individual storm 
intensity is changing through time. Any increase in cumulative 
yearly storm intensity and potential destructiveness, therefore, is 
due to the increasing number of storms and not due to any increase 
in the intensity of individual storms. 



I. Introduction 

It is important to be able to statistically characterize the distribu- 
tion of the number of hurricanes in the North Atlantic, especially if this 
number is increasing or the intensity of hurricanes is changing. Sev- 
eral studi es hay e examined this. The most important recent paper is 



Emanuel (120051 ). in which he argued that hurricanes in the North At- 
lantic have become more destructive over the past 30 years. To measure 
potential "destructiveness", he developed a measure called the power 
dissipiation index, which is a function of the cubed wind speed of a 
storm over its lifetime (see Sec. [2] for a precise definition). In his 
original paper, this index was not just a measure of a single storm's 
intensit y, bu t a cumulative index over all the storms during the year. 



Peilke (120051 ) and Landsea (120051 ) criticized the data analysis method 
used to demonstrate that the index was increasing; by pointing out that 
the smoothing method used on the raw time series data was slightly 
flawed, that errors in the observations should lead to a less certain 
statement about increases, and that the wind speed adjustments used 
by Emanuel were too aggressive. 

Our observation is that the lumping together of all the storms within 
a year has lead to a different interpretation of what exactly is increas- 
ing: storm number or (a function of) windspeed. Other explanations of 
Emanuel's findings may be that the number of cyclones has remained 
(distributionally) constant, but that average storm intensity has in- 
creased. Or it may also be that the number of cyclones has increased 
but that the intensity on individual storms has remained constant, 
or even decreased. Other combinations are, of course, possible: both 
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storm frequency and individual storm intensity might have increased. 
We examine these scenarios below. 

A fi rst step in such an analysis was taken by Eisner and Bossak 
( 120011 ) . who examined the climatology and seasonal modeling of hur- 
ricane rates using a series of Bayesian hierarchical models. Using this 
modern approach allows us to easily make probability statements about 
important parameters of storm correlates and t o specify the uncertainty 
of future predictions. Eisner and Jagger ( 20041 ) continued the Bayesian 
modelling line by controlling, in their models, for the influence of the 
cold tounge index and the North Atlantic oscillation index. They found 
that both of these indices were well correlated with the mean hurricane 
numb er. We also use these indicies in our models below. Eisner et al. 
( 120011 ) investigate th e rela tionship between ENSO and hurricane num- 
bers. Hoyos et al. (l2006h examine these and other factors that may 
contribute to in crease s in the mean frequency of hur ricanes. 

Eisner et al. ( 2004 ) and Jewson and Penzer ( 20061 ) examine whether 
there were shifts, or change points, in the statistical distribution of hur- 
ricane numbers. Both groups of authors did find likely changes, namely 
around 1900, the mid 1940s, mid 1960s, and the mid 1990s. These shifts 
may have been due to actual changes in physical mechanisms (such as 
large-scale shifts in the atmospheric or oceanic circulations) or they 
may be due to changes in measurements, though all authors agree that 
the changes are probably a combination of both. We take this topic up 
below, but we do not seek to answer why these changes take place, or 
even if they are certain. It is clear enough, however, that the data has 
changed in character through time. Thus, we build our models using 
different ranges of data in an attempt to incorporate this uncertainty. 

Our approach is different from previous analyses in two ways: (1) 
we hierarchically model the number of tropical storms and then the 
chance that hurricanes arise from them, as opposed to directly model- 
ing t he number of hurrican es (or hurricane land falls); and (2) in line 



with Webster et al.l (120051 ) we characterize the distribution of storm 
intensity within a given year and ask whether this distribution changes 
through time. We investigate storm intensity by measuring storm life- 
time in days, storm track length, and Emanuel's power dissiptation 
index applied to invididual storms. 

We use the hurricane reanalysis datab ase (HURDAT) Jarvinen at 
al. ( 20021 ) and, e.g. Landsea et al. ( 2004 ). This database contains six- 
hourly maximum sustained 1-minute winds at 10 m, central pressures, 
and position to the nearest 0.1° latitude and longitude from all known 
tropical storms from 1851-2006. A cyclone was classified as a "hurri- 
cane" if, at any time during its lifetime, the maximum windspeed ever 
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met or exceeded 65 knots. Obviously, this cutoff, though historical, is 
somewhat aribitrary and other numbers can be used: we discuss this in 
greater detail below. To investigate the realtionship of North Atlantic 
tropical storms with ENS O, we al so use the cold tongue index (CTI) 
Deser and Wallace! (119901 ) and e.g. ICail (120031). And we u se the North 
Atlantic oscillation index (NAOI) from I Jones et al.l (119971 ). 

Section 2 lays out the statistical models and methods that we use, 
Section 3 contains the main results, and Section 4 presents some ideas 
for future research. 



2. Methods 

We adopt, as have many before, Bayesian statistical models. An 
important advantage to these models is that we can make direct prob- 
ability statetments about the results. We are also able to create more 
complicated and realistic models and solve them using the same numer- 
ical strategy; namely, Gibbs sampling. We do not go into depth about 
the particular methods involved in forming or solving these models, as 
readers are likely familiar with these metho ds nowadays. There are 



also many excellent references available, e.g. iGelman et al.l ( 120031 ). 

It is important to control for factors that are known to be related, 
or could cause changes in, the frequency of tropical cyclones and storm 
intensity. There are many such possible variables, but we choose, for 
ease of com pariso n, the same indicies as cited in the paper by Eisner 
and Jagger (120041 ) . These factors are the cold tongue index and the 
North Atlantic oscillation index. Readers are encouraged to refer to the 
original sources and references therein to learn about these indicies. 

2.1. Number of storms. Most statistical analysis focuses on the num- 

ber of h urricanes or the sub s et of l andfalling hurricanes, e.g. lElsner and Jagger 
( 2004 ); Eisner and Bossakl ( 200ll ). The approach here is different. We 
first model the number of tropical cyclones and then model whether or 
not, for any given tropical cyclone, a hurricane evolves from it. Specifi- 
cally, we do not separately model the frequency of both hurricanes and 
cyclones, as doing this ignores the relationship of how cyclones develop 
into hurricanes. 

We suppose, in year i of n years, that the number of storms is well 
approximated by a Poisson distribution as in 

(1) Si ~ Poisson(Aj) 

where Aj describes the mean (and variance) of the number of storms. 
It is of primary interest to discover whether this parameter is changing 
(possibly increasing) through time, controlling for known important 
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meterological and oceaniographic variables. Eisner and Jagger (120041 ) 
developed this same model for the number of hurricanes (and not cy- 
clones per se). Here we adapt it to the number of cyclones, and add in 
the possibility that the parameter A changes in a linear fashion in time 
(we also, for ease of reference, adopt Eisner and Jagger's notation). 
Thus, we further model Aj as a function of the CTI and the NAOI, and 
allow the possibility that Aj changes linearly through time. We use the 
generalized linear model 

(2) log(Ai) =p + Pit + /9 2 CTI 4 + /3 3 NAOL, + &CTI, x NAOI, 
The prior for (3 is 

(3) /3 ~ MVNipX 1 ) 

that is, (3 = (/?o,/3i, ■ ■ ■ , Pa) is distributed as multivariate Normal, and 
t = 1, . . . ,n. We take noninformative values on the hyperparameter 
a = (0, 0, . . . , 0). Two priors were tried for b" 1 = 0, which is equivalent 
to the standard improper flat prior, and a general covariance matrix 
where b was given values on the off diagonals to allow for correlations 
between the betas; both priors gave nearly equivalent results both in 
this model and those below. We report results on the standard prior 
(technically, b~ l is the precision and b^ 1 = is the default "ignorance" 
prior in the software, which is described in Section 3). Thus, if the 
posterior, for example, Pr(/?i > 0|data) is large then we would have 
confidence that the mean number of storms is increasing. 

We will also be interested in two other posterior probability distribu- 
tions, useful for model checking: p(Aj|data) = Jf° p(Aj|/3, data)p(/3|data)<i/3, 
which is the posterior of the mean number of storms in year i integrat- 
ing over the uncertainty in (the multidimensional) Pi, and p(gi\ data) = 
/o°°P(^i|^i! data)p(Ai|data)c/Aj, where the later represents the posterior 
predictive density (a guess) of the number of storms #j in year i inte- 
grating out the uncertainty in Aj represented by p(Aj|data): to be clear, 
p((?i|data) is the distribution we would use if we wanted to make a fore- 
cast for Si given t, CTIi, and NAOIi, integrating over the uncertainty 
we have in P and Aj. These two distributions are not analytic and in 
practice we sample from them using Monte Carlo methods, which are 
described below. 

Once a tropical storm develops it, of course, has a chance to grow 
into a hurricane. If there are s tropical cyclones in a year the number 
of hurricanes is constrained to be between and s. Thus, a reasonable 
model for the number of hurricanes hi in year % given Sj is 

(4) hi\si ~ Binomial(sj, #j) 
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Now, if we assume that the prior on B{ is 

(5) 6i ~ Beta(a, b) 

with a = b = 1 (which represents a flat prior) then the posterior of 0, 
is 

(6) Oi\hi, Si, a, b ~ Beta(/ij — 1, s, — /ij — 1) 

which makes calculations easy. This same model could of course be 
adapted to define hurricanes as "major storms" using cutoffs greater 
than 65 knots, say greater than 114 knots for category 4 storms. We 
could then estimate how many tropical cyclones evolve to these more 
destructive storms (we plan to address this more completely in future 
work, though we make a few comments below). 

It is possible, however, as with Aj, that 9i is dependent on CTI and 
NAOI and that it changes through time. To investigate this, we adopt 
the following logistic regression model 



(7) log [y^J.) = h + ^ + ^ CTli + &NAOL, + /3 4 CTIi x NAOI, 
where we again let 

(8) P~MVN(a,b) 

where a, b are as before. And again, if, for example, Pr(/?i > 0|data) 
is large then we would have confidence that the number of storms 
that turn into hurricanes is increasing. In practice, however, we first 
transform the posteriors by exponentiation, thus the posteriors have 
the natural interpretation of odds ratios for the increase (or decrease) 
in theta. So we would instead look at Pr(exp(/?i) > l|data), and if this 
is large, we would have good evidence that the mean rate of cyclones 
converting to hurricanes is increasing, i.e the odds of converting are 
greater than 1. 

2.2. Measures of intensity. It may be that the frequency of storms 
and hurricanes remains unchanged through time, but that other char- 
acteristics of these storms have changed. One important characteristic 
is intensit y We define three m easures of intensity, in line with those 
defined in Webster et al. fl2005h : (1) the length m, in days, that a storm 



lives; (2) the length of the track (km) of the storm over its lifetime; 
and (3) the power dissipation index as derived by Emmanuel, though 
here we apply this to each cyclone individually. 

m was directly available from the HURDAT reanalysis: we approx- 
imate the number of days to the nearest six-hours. Track length was 
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estimated by computing the great circle distance between succesive six- 
hour observations of the cyclone, and summing these over the storm 
lifetime. The power dissipation index (PDI) is defined by 

(9) PDI = fv^Jt 

Jo 

where V^ ax is the maximum sustained wind speed at 10m, and T repre- 
sents the total time that the storm lived. Practically, we approximate 
the PDI — up to a constant — by summing the values (Knax/100) 3 at 
each six-hour observation. The PDI is a crude measure of the strength 
of the potential de struct iveness of a tropical storm or hurricane, as 



cited by Emanuel (120051 ). Other than this measure, we say nothing 
directly about storm destructiveness (in terms of money etc.). 

It was found that log transforms of these variables made them much 
more managable in terms of statistical analysis. Transforming them led 
to all giving reasonable approximations of normal distributions; thus, 
standard methods are readily available. 

For all three of these measures, we adopt a hierarchical modelling 
approach because we are interested if the distribution within a year 
of these measures changes through time. Let x^j be any of the three 
measures (m, track length, or PDI) for storm % in year j. Then we 
suppose that 

(10) log a;*,,- ~iV(7i,i,x) 
where 

(11) 7ij - = fa + + AiCTI, + AjNAOI, + faCTIj x NAOL, 
and, as before, where 

(12) /?~MWV(/i,r) 

and fi = 0, and x, T -, which are assumed independent of fi, have inverse 
Gamma(0.005, 0.005) flat priors. This model is, of course, similar to 
those above, except that we seek to characterize the distribution of 
intensity within each year and see how that distribution mean might 
change in a (linear) way. 



3. Results 



AH computations were carried out in the R statisti c al sys tem lR Development Core Team 
(j2005l ) using the MCMC package [Martin and Quinnl (j2007l ) on a Fedora 
Core 6 platform. Models were fitted using Gibbs sampling. The first 
1000 simulations were considered "burn in" and were removed from the 
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analysis: 10,000 additional samples were calculated after this and used 
to approximate the posterior distirbutions. 

Data for all the measures we use was available from 1851, but we 
use only data from 1900 onwards. The data from before this date, 
as is well known, is suspect enough to cast suspicion on any results 
based on them. It is also not clear that a strict linear model over 
the entire period of 1900-2006 would best fit these data a s observation 



and in strument changes through that time have changed lElsner et al. 



(120041 ) . So we adopt the practice of computing each model over three 
different time periods: once for the entire period 1900-2006; the second 
for dates between 1950-2006; and the third between 1980-2006. These 
choices are somewha t arbi tr ary, but in line with the c hange-point re- 
sults of Eisner et al. ( 2004 ): Jewson and Penzer ( 20061 ). Other choices 
are easily made, however, and we have found that our results are ro- 
bust to changes in these exact start times. This approach also lets us 
check whether a linear model for increase/decrease of the parameters 
through time is reasonable. We do not investigate more complicated 
models, such as linear change point regression models, here. Lastly , 
we do not adju s t the obs e rved wind speed in any way IPielkd fl2005h : 



Landseal (120051 ): lEmanuell (120051 ). 



3.1. Number of storms. The top two panels of Fig. [T] shows the 
time series plots of s and h. There does appear, to the eye, to be an 
increase in s in the past two decades and perhaps a smaller increase in 
h, in line with what others have found. 

Figure [2] shows the posterior distributions from the model ([2]). Table 
1 gives the summary statistics for this model. In each case, and in 
all future figures, the solid line represents the model using all data 
from 1900-2006; the dashed line represents the model using data from 
1950-2006; and the dotted line represents the model using data from 
1980-2006. Regardless of the data used, there is good evidence that 
Pr(/9i > 1 data) ~ 1, which implies that the mean number of storms 
has increased through time. However, using only the most recent (1980- 
2006) data does not give a very certain estimate for /3\, which can be 
see by noting that the probability for (3i taking any particular value is 
distributed over a large range of possible values. 

The parameter is the "intercept" in model ([2]), and is the estimate 
of Xi when the other parameters, time included, are set to 0. Thus its 
particular value does not directly mean much. We can note that the 
intercept for A« has shifted to higher values when considering the two 
later data sets, though we can not rule out that there has actually been 
change because the three distributions show great overlap. 
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Figure 1. The number of storms s and hurricanes h, 
and the ratio h/s, for the North Atlantic from 1900-2006. 



There is also strong evidence, as Eisner and Jagger found, that the 
CTI (P2) is important in estimating Xf. greater CTIs lead to smaller 
As, and therefore to a smaller probability that the mean number of 
storms will be high; or more plainly, greater CTI means fewer storms. 
It appears that this realtionship has strengthened in later years (1980- 
2006), as the distribution has shifted to smaller numbers. 
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Results for the NAOI ((3$) are more mixed. Data from 1900-2006 
and 1950-2006 do not give high probability that NAOI is influential. 
But when only the 1980-2006 data is used, the probability is high that 
when NAOI is greater that A will be smaller, and again, a smaller 
probability that the number of storms will be high. The same goes for 
the interaction between CTI and NAOI (fa), though even using the 
1980-2006 data does not give conclusive results. 

The Table presents the same data as the figures, but only for the 
1980-2006, in tabular form, so that readers can read a best (assuming 
absolute error loss) estimate of each (3 at the median (50%-tile). 95% 
credible intervals can also be read from the table: the values from 
the 2.5%- to the 97.5%-tiles. The last column gives the estimated 
probability that P(/3 > 0|data). For fli, if this number is large, it 
means we have great confindence that increases are taking place. For 
the other (3s, if this number is near 0, it means that we are confident 
that P((3 < 0|data) (= 1 - P{j3 > 0|data)). 

These, and the following results, are only as good as the model. To 
check model quality, we present Figs. [3J and HI plotted for the 1980- 
2006 data. Fig. [3] is a plot of the two posterior distributions p(Aj|data) 
and p(#j|data). Overplotted on both of these distributions is the actual 
value (in black dots) of Sj. Recall, that p(Aj|data) is our best guess of 
the mean number of storms, and not the number of storms themselves. 
This is the right-hand side of the figure: p(<7j|data). To calculate these, 
we use the same MCMC results that came from the generation of the 
posteriors of j3, and plug these into (|2J) (along the with the appropriate 
values of t, CTI iy and NAOIj) and then solve for Aj. We first sort the 
actual values of s and plot density estimates of the posteriors, adding 
an arbitrary number to each so that they may all fit on the same 
plot. Fatter plots indicate more uncertainty as to the exact value of Aj; 
narrower plots indicate the opposite. Most of the dots are somewhat 
near the peaks of these posteriors, but certainly not all. 

The right-hand side shows the posterior of the best guess gi. These 
posteriors are got from, at each of the 10,000 samples of each Aj, simu- 
lating 50 Poisson random variables with mean A^: then, the 50 x 10, 000 
values are used in building density estimates of the posterior predicitve 
gfj. Notice that these plots are jagged, which is due to the discrete 
nature of Qi and Sj. Here, to give some relief to the figure we estimated 
the density and took the cube root of each frequency value; which has 
the effect of exaggerating the high (mode) values so they are easier to 
see. These posteriors are certainly more spread out, as expected. For 
example, if g ~ Poisson(12) then a 90% prediction interval is 8 to 17. 
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Table 1. Common quantiles of the model parameters 
and the posterior probability that these parameters are 
greater than for 1980-2006 data only. The exception is 
model (EP where the probability is that the parameter is 
greater than 1 (because these are odds ratios). 



Parameter 


2.5% 


50% 


97.5% 


P(P > 0|data) 


s 


Pi 


0.011 


0.026 


0.041 


0.9996 


CTI 


-0.61 


-0.40 


-0.19 


< le-5 


NAOI 


-0.40 


-0.22 


-0.059 


0.007 


CTIxNAOI 


-0.40 


-0.16 


0.077 


0.095 






h/s 




(PQ3 > l|data)) 


Pi 


0.997 


0.9999 


1.002 


0.47 


CTI 


0.46 


0.71 


1.098 


0.066 


NAOI 


0.61 


0.90 


1.29 


0.30 


CTIxNAOI 


0.68 


1.14 


1.90 


0.69 






Cat IV/s 




(P(P > l|data)) 


Pi 


1.0 


1.004 


1.009 


0.99 


CTI 


0.13 


0.30 


.64 


0.0001 


NAOI 


0.39 


0.84 


1.69 


0.33 


CTIxNAOI 


0.19 


0.48 


1.20 


0.06 






log(m) 






Pi 


-0.006 


0.003 


0.011 


0.74 


CTI 


-0.20 


-0.075 


0.053 


0.13 


NAOI 


-0.13 


-0.023 


0.085 


0.34 


CTIxNAOI 


-0.07 


0.080 


0.23 


0.86 




log(track leng 


th) 




Pi 


-0.0073 


0.004 


0.012 


0.76 


CTI 


-0.25 


-0.079 


0.091 


0.18 


NAOI 


-0.22 


-0.081 


0.063 


0.13 


CTIxNAOI 


-0.18 


0.019 


0.22 


0.57 






log(PDI) 






Pi 


-0.012 


0.0021 


0.022 


0.58 


CTI 


-0.65 


-0.36 


-0.056 


0.01 


NAOI 


-0.34 


-0.09 


0.16 


0.24 


CTIxNAOI 


-0.31 


0.043 


0.40 


0.59 



However, it is clearer from this picture that the model is not behaving 
too badly. 

More practically, a forecaster would pick a single number g^ as a 
guess, plus a measure of uncertainty of this guess. Fig. @] shows the 
actual values Sj (slightly jittered to separate the points) by the median 
of gfj|data. A one-to-one line is also shown. Thick black lines are the 
50% credible interval, and the dotted thin lines are the 90% credible 
interval. All but 7 of the points are within the 50% prediction intervals, 
and all but one (2005) are within the 90% prediction intervals. This 
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gives us confidence that the model we have chosen is representative of 
the real ocean-atmosphere; though it is, of course, not perfect. 

The bottom panel of Fig. [I] shows the ratio hi/ s-i which does not 
give much indication of changing through time. Under this assumption, 
and using the well-known properties of the posterior (jSJ) for the data 
1900-2006, we estimate that the mean fraction of converting storms is 
^ = = §§t = 0-^6, and standard deviation of 0.01 (calculation 

not shown). That is, once a tropical cyclone forms, there is a 56% 
chance that it will evolve into a hurricane. Results for 1950-2006 are 
58%, and for 1980-2006 are 55% (with similar standard deviations). 

But applying model ([6]) gives Fig. [51 summarized also in Table 1. 
The rate at which storms become hurricanes does not appear to change 
through time, evidenced by the probability of f3\ not near 1 is small: 
recall that these results are in the form of odds ratios, except for /? , 
which is in the form of odds. Actually, there is some evidence, using 
the data from 1950-2006, that mean rate at which storms evolved was 
actually less, though the effect, if real, was quite small. The best 
estimate, using 1980-2006, is that the odds of cyclone evolving into a 
hurricane are 1 per year; which is, of course, no change at all. There is 
weak evidence that CTI contributes negatively, in the sense that when 
CTI increases, the odds of a cyclone becoming a hurricane decreases, 
by the odds of about 0.7 times per CTI unit. NAOI does not appear 
important (most of the probability mass is around the odds ratio of 1). 
The interaction of CTI and NAOI does not appear important. 

The conclusion is that there is good evidence that the number of 
tropical cyclones has increased, but that the chance that these cyclones 
become hurricanes has remained constant, even after controlling for 
CTI and NAOI. This also means, of course, that it is likely that the 
number of hurricanes has also increased. 

We repeated this analysis (plots not shown, though results are in 
the Table) for category 4 and above hurricanes using the 1980-2006 
data (defined as a cyclone which had winds meeting or exceeding 114 
knots at any time during its lifetime), and find that the odds that 
the rate these storms evolve from tropical cyclones has increased by 
the odds of only 1.004 times per year (95% credible interval 1.000 to 
1.009). Much of this increase is due to 1999, when 42% of the cyclones 
evolved into major storms, whereas the average rate over the remaining 
years is about 10%. Now, it is well known that standard regression 
models, such as this one, are quite sensitive to these sort of "outliers" , 
so this result should be taken with a very large grain of salt. Too, the 
effect disappears when the data from 1950-2006 is used (95% credible 
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interval 0.998 to 1.002). There is good evidence that CTI contributes 
negatively, in the sense that when CTI increases, the odds of a cyclone 
becoming a major storm decreases, by the odds of about 0.3 times per 
CTI unit. NAOI does not appear important (most of the probability 
mass is around the odds ratio of 1). The interaction of CTI and NAOI 
does appear important, though this is in large part driven by the CTI 
values. 

3.2. Measures of intensity. Figure [6] shows the time series boxplots, 
for each year, of log(m), log(track length), and log(PDI). The boxplot 
gives an indication of the distribution of each measure within each year. 
There is no apparent trend in the sense that, say, the medians show no 
systematic direction, and the other quantiles appear distributed around 
a central point. If this graphical view holds under modeling, it means 
that cyclone distribution of intensity has not changed through time. 

We now apply model (fTTl) to each of these measures. Figs. 019] and 
Table 1 summarize the results. 

For none of the measures, using any of the data time sets, does there 
appear to be substantial evidence that the mean of the distribution of 
these measures changed though time: the posteriors for j3i in each case 
have most of their mass around 0. Except perhaps for (log) PDI using 
the 1950-2006 data set. Here there is good evidence that mean (log) 
PDI has decreased over this time period because most of the probability 
mass is at values less than 0. (P(j3 < 0|data) = 0.998). Recall that 
this is the distribution each individual storm's PDI within a year, and 
not a yearly (additive) summary. 

For the (log) number of days or (log) track length, there is not much 
evidence that CTI or NAOI are influential, evidenced by the high prob- 
abilities of values around 0. This is not so for (log) PDI, where there 
is some evidence that higher CTI and (and maybe higher NAOI) lead 
to smaller mean (log) PDIs. The interaction between CTI and NAOI 
does not appear important. 

4. Conclusions 

We find that there is good evidence that the number of tropical 
cyclones in the North Atlantic basin have increased in the past two- 
and-a-half decades. This result stands even after controlling for CTI 
and NAOI. These results are of course conditional on the model we 
used being adequate or at least it being a reasonable approximation to 
the data. Diagnostics plots (Figs. [3] and HJ) indicate that the model 
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performs well, though of course not perfectly, at predicting the histori- 
cal data. We make no predictions about future increases as it would be 
foolish to extrapolate the simple linear model we used into the future. 

We also found that the rate at which tropical cyclones become hur- 
ricanes does not appear to be changing through time, nor is it much 
influenced by CTI or NAOI. There is weak evidence that the mean 
rate at which major (category 4 or above) storms evolve has increased 
through time, though the increase is very small, and it is just as likely 
the model used to assess this rate is inadequate. 

We find almost no evidence that the distribution of individual storm 
intensity, measured by storm days, track length, or individual storm 
PDI, has changed (increased or decreased) through time. Any increase 
in storm intensity a t the conglomerate yearly level, as for example found 
by lEmanuell (120051 ). is likely due to the increased number of storms and 
not by the increased intensity of individual storms. We also repeated 
our analysis on the distribution of each storm's (log) maximum wind 
speed over its lifetime and came to the same conclusion as with the 
other measures of intensity. 

Much more exact work can be done. A model similar to that above 
should certainly be used for storms across all ocean basins for whic h 
data is available, as was recently done by Webster et al. (120051 ). 
Too, more sophisticated models could be used. For example, spa- 
tial Bayes ian models such as those developed by Wikle and Ander- 
son (120031 ) for estimating tornado frequency change could be used for 
tropical cyclones. This is not an easy task because tornadoes, in that 
model, were treated as point objects, and, of course, hurricanes vary 
in intensity over vast spatial lengths. The statistical characteristics of 
individual tropical cyclones could be better addressed, by asking how 
the change in intensity (by the three measures given above, and by 
others such as pressure), changes through storm lifetime. And a more 
complete model state-space like model for the multivariate measure of 
intensity (days, track length, PDI) that takes into account the correla- 
tion between these dimensions could certainly cast more light on how 
individual storm intensity may or may not have changed through time. 
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Figure 2. The posterior distributions for the parame- 
ters in model ([2]). In each case, and in all future figures, 
the solid line represents the model using all data from 
1900-2006; the dashed line represents the model using 
data from 1950-2006; and the dotted line represents the 
model using data from 1980-2006. Regardless of the data 
used, there is good evidence that Pr(/?i > 0|data) ~ 1, 
which implies that the (rate of the) number of storms is 
increasing through time. The other parameters are dis- 
cussed in the main text. Posteriors with most of their 
mass around have little to no effect on the outcome. 
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Figure 3. The (sorted by number of storms) posteriors 
p(Aj|data) and p(#j|data) for the tropical cyclones during 
1980-2006. Black dots show the actual values. See the 



text for details on how these plots were constructed. 
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Figure 4. The medians (predictions) of each p^data 
with 50% credible intervals (solid lines) and 90% credible 
intervals (thin, dotted lines), plotted by the actual values 
of Si. Each Si was slightly jittered to aid the eye. 
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Figure 5. As in Fig. (j2D except for model (El). The 
results here differ in that the posteriors are in terms of 
odds ratios for all but the intercept, which is in odds. 
Thus, posteriors with most of their mass around 1 have 
little to no effect on the outcome. 
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Figure 6. The time series of boxplots, for each year, 
of log(m), log(track length), and log(PDI). There is no 
apparent trend. 
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Figure 7. As in Fig. flSJ) except for model ([II]) for the 
logged number of storm days. 
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Figure 8. As in Fig. fl2J) except for model (p]) for the 
logged track length. 
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Figure 9. As in Fig. fl2J) except for model (p]) for the 
logged PDI. 



